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Clustering of matter in waves and currents 
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The growth rate of small-scale density inhomogeneities (the entropy production rate) is given by 
the sum of the Lyapunov exponents in a random flow. We derive an analytic formula for the rate in 
a flow of weakly interacting waves and show that in most cases it is zero up to the fourth order in 
the wave amplitude. We then derive an analytic formula for the rate in a flow of potential waves and 
solenoidal currents. Estimates of the rate and the fractal dimension of the density distribution show 
that the interplay between waves and currents is a realistic mechanism for providing patchiness of 
pollutant distribution on the ocean surface. 

PACS numbers: 47.27.Qb, 05.40.-a 



Random compressible flows produce very inhomoge- 
neous distribution of density, see fil ^_ E 0, IE IE 0, [EH] 
for theory and [IE EH, EE EE EE EE El | for experiments. 
Here we study the growth of density inhomogeneities at 
small scales, where the flow can be considered spatially 
smooth. It can then be characterized by the Lyapunov 
exponents whose sum is the logarithmic rate of change of 
an infinitesimal volume element, that is minus the den- 
sity rate of change A. It is called also the entropy pro- 
duction rate or the clustering rate. Since contracting 
regions contain more statistical weight than expanding 
ones, A is generally positive in a random compressible 
flow @, IE El HI] (an analog of the second law of thermo- 
dynamics). As a result, density grows on most trajecto- 
ries. In the limit of infinite time, density concentrates on 
a constantly evolving fractal set characterized by a sin- 
gular (Sinai- Ruelle-Bowen) measure [J, EE EE EH ■ Our 
goal here is to establish what determines the rate A in 
fluid flows with waves, particularly, on liquid surfaces. 
Patchiness in the distribution of litter on the surface of 
lakes and pools and of oil slicks and seaweeds on the sea 
surface is well-known empirically while there is no the- 
ory that describes it. We think that in many situations 
patchiness is a signature of a fractal measure forming on 
the surface. Our purpose here is to estimate how fast this 
fractal set is formed and what is its fractal dimension. 

Any velocity field can be separated into solenoidal (in- 
compressible) and potential components. Surface flows 
can be compressible, even for incompressible fluids. For 
example, underwater turbulence produces compressible 
surface currents that lead to fractal distributions of sur- 
face density [Tl|, EE EE EE El) HE] • However, underwater 
turbulence is relatively rare in natural environment (be- 
cause of stable stratification) and large-scale currents are 
usually incompressible. Compressible component of the 
surface flows is then provided by potential waves. Linear 
waves just oscillate, net effects are produced by nonlin- 
earity. Every running plane wave provides for a (Stokes) 
drift proportional to the square of the wave amplitude. A 
set of random waves provides for the Lyapunov exponents 
proportional to the fourth power of the wave amplitudes 
yet the sum of the exponents —A is found to be zero for 



purely potential waves with Gaussian statistics (nonzero 
rate appears only in the sixth order in wave amplitudes, 
i.e. it is so small as to be practically unobservable in most 
cases) [I]. Here we use the recently derived general for- 
mula for the entropy production rate (2l| and show that 
the account of wave interaction (which makes the statis- 
tics weakly non-Gaussian) does not bring non-zero A in 
the fourth order in wave amplitudes. We then suggest 
that in many situations (particularly on liquid surfaces) 
the growth of density inhomogeneities is due to an in- 
terplay between potential waves and solenoidal currents. 
For such flows, we calculate A and the fractal dimension 
of the resulting measure and consider different limits. 

In the velocity field v(t, x), the trajectory X(t, x) sat- 
isfies the equation <9 t X = v(t, X) with the initial condi- 
tion X(0,x) = x. The rate of density change along the 
trajectory averaged over x is given by plj : 



A = - lim (w(t,X)) 



dt(w(Q,x)w(t,X)) , (1) 



with w = V • v. This is a generalization of the Kawasa- 
ki formula [22| (obtained in the context of statistical 
physics) to time-dependent flows with a steady statistics. 

For a general flow, it is impossible to relate the La- 
grangian integral |T]) to the velocity spectra or correlation 
functions given usually in the Eulerian frame. However, 
for low-amplitude waves, fluid particles shift little dur- 
ing a period, which allows for an analytical treatment. 
Considering waves with the dispersion relation and 
packets with both the wavenumber and the width of or- 
der k, we estimate the correlation time of w as 1 and 
the correlation length as k~ x . The deviation X(t, x) — x 
during t ~ fi^ 1 is e = few/fife *C 1 times smaller than 
k^ 1 . We now expand |T]) near x up to e 4 : 
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dhvP(t 2 )) (2) 
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All quantities here are taken at the same point in space. 

We start the consideration of @ from the simplest case 
when the flow is solely due to weakly nonlinear waves. 
The normal coordinates of such waves satisfy the equa- 
tion dtdk = —iSH./5a k while the velocity Fourier com- 
ponent is assumed to be given by Vk = Ak(ak — <iL k )- 
Here H. is the wave Hamiltonian, which can be expanded 
in wave amplitudes as follows [23|]: H = J dkfl k \ak\ 2 + 
\ J dkus (Vi23ai(i2 a 3 + c.c.) + • ■ ■ . We do not write ex- 
plicitly here other (third and fourth-order) terms since 
they will not contribute A up to ~ e 4 . We use throughout 
the shorthand notations V123 — Vi23<5(fci — &2 — ^3) and 
V 123 = V(k u k 2 ,k 3 ), Q(±ki) = n ±i and A(ki) = A,. 

One derives the clustering rate up to e 4 using a stan- 
dard perturbation theory for weakly interacting waves 
[23 | . The first term in @ is the time integral (the zero- 
frequency value) of the second moment. At the order 
e 2 , the second moment in the frequency representation is 
proportional to the delta function: (a*(k, u>)a(k', u>')) — 
{2i:) d+2 n(k)8(uj-n k )5(k-k')8{uj~uj') . A finite width 
over u> and a finite value at u) — appear either due to fi- 
nite linear attenuation (the case considered in [24[ ) or due 
to nonlinearity in the second order of perturbation the- 
ory (which gives e 4 and is considered here). The second 
term in ([2]) is the triple moment which appears in the first 
order of the perturbation theory and the last two terms 
contain the fourth moment which is to be taken at the 
zcroth order (i.e. as a product of two second moments). 
Straightforward calculations then give for weakly nonlin- 
ear waves the e 4 contribution: 



A = Re 



dk2<ik 3 
(27r) 2d 



S(Q 2 - n 3 )n(fc2)n(fc 3 ) 



dki 

(2n) a 



V ™ Vs - 12 ^^(27r) 3d+1 \A l .k 1 \ 2 ^- (3) 



(27r) 2d 

^2 



^ I (A 3 • k 3 )(A* 2 ■ k 3 ) - (A 3 • fe 2 )(A| • fe) 



{A\ •fc 1 )(A 2 -fc 2 )(A*-(fc2 + fci)) 



(4) 



(5) 



The common factor S(fl2 — fl3)n(k2)n(k3) tells that we 
have here the contribution of two pairs of waves with 
the same frequencies. All three terms are generally 
nonzero (and positive) when the dispersion law is non- 
monotonic or non-isotropic so that Q2 = ^3 does not re- 
quire fc 2 = &3- In most interesting cases, however, flk is 
a monotonous function of the modulus k so that k% = fc 3 - 
Let us show first that wave interaction does not con- 
tribute A in this case. Indeed, the first two terms, (13l4p . 
that came out of the first two terms of ([2]), are propor- 
tional to the difference, V^* 13 — V3_i 2 , between the ampli- 
tude of decay into a wave with k\ and confluence with a 
wave with — k\. Interaction coefficients for fc 2 = &3 have 
rotational symmetry and are thus functions of wavcnum- 
bers so that V213 - Va_i2 = V213 - V312 = V212 - V212 = 0. 

The last term (JSJ comes from the last two terms of © 
and does not contain the interaction coefficient V. This 



term is due to nonlinear relation between Eulerian and 
Lagrangian variables rather than due to wave interaction. 
We can compare JS]) with the growth rate of the squared 
density for non-interacting waves, see (12) in [8] writ- 
ten there in terms of the energy spectrum, E a ^{k 1 uo) — 
2TtA a {k)A*P{k)[n{k)5{u)-tt k )+n{-k)8{Lu + ri- k )]. The 
comparison shows this part of our logarithmic growth 
rate being exactly half the growth rate for the second 
moment as it should be for a short-correlated flow . In- 
deed, the process of creation of density inhomogeneities is 
effectively short-correlated since the time it takes (1/ fi^e 4 
or longer) exceeds the correlation time of velocity diver- 
gence in the Lagrangian frame, 1/flk- For monotonous 
fi(fc), ([5]) is nonzero only if the polarization vector A k 
is neither parallel nor perpendicular to k i.e. contains 
both solenoidal and potential components. This is not 
the case for most waves in continuous media. We thus 
conclude that for most common situations (in particular, 
for sound or surface waves) the entropy production rate 
is zero in the order e 4 . Note that for surface waves, the 
canonical variables are elevation r](r,t) and the potential 
</>(r, z = rj,t) which are related to the surface velocity by 
a nonlinear relation v = V0(r, 77, t). Expanding it in the 
powers of 77, one can show that this extra nonlinearity 
does not contribute A in the order e 4 ;25||. We find it re- 
markable that the flow of random potential waves is only 
weakly compressible (i.e. the senior Lyapunov exponent 
is much larger than the sum of the exponents). 

Therefore, we consider now the clustering rate in the 
presence of solenoidal currents and potential waves, the 
situation most relevant for oceanological applications 
[2^ |. Consider the solenoidal flow u weakly perturbed 
by potential waves with v. To derive A in the lowest 
(second) order in e = kv/Qk, we neglect the contribution 
of v into X in (fT]) and assume d t X(t,x) w u(t,X(t,x)). 
In this order, w = V-t? is Gaussian and one may integrate 
by parts: (w(Q, x)w(t, X(t, x))) = J dt'dx'$(t' ,x' - 
x)(5w(t,X(t,x))/5w(t',x')}. Here $(t' - t,x' - x) = 
{w(t, x)w(t', x')) is the Eulerian correlation function and 



/>oo 

/ dt(<&[t,J{t)}), J(t) = X(t,x)-x. 
Jo 



(6) 



Waves and currents are considered statistically indepen- 
dent in this order. Using the spectrum, k a k$E^f = 
k 2 E k , we can express ®(t, r) = {2n)~ d J k 2 Ek cos(fc • r — 
flkt)dk and rewrite ^ as a weighted spectral integral: 



A = (2tt 



,-d 



k 2 Ekn{k)dk, 



fi(k) = / (cos [k ■ J(t) - n h t]) dt. 
Jo 



(7) 
(8) 



The spectral weight fi(k) is the Lagrangian correla- 
tion time of the fc-harmonic of w and is expressed via 
the characteristic function of the particle drift J it). 
Without currents, (|7l8p reproduce the first term of ([2]) 
since only the zero- frequency wave contributes. Al- 
ready a steady uniform current u contributes the clus- 
tering rate in the order e 2 if there are waves whose 
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Doppler-shifted frequency is zero in the current refer- 
ence frame: A = (27r)~ d J k 2 E k S(Q k — k ■ u)dk. Sim- 
ilar Cherenkov resonance has been noticed before for 
diffusivity 27j. In the rest of the paper we consider 



the fluctuating part of current velocity characterized by 
the rms velocity u 2 = (u 2 ) and the correlation time 
t = J(u a (0, x)u a (t, X(t, x))dt/v,Q. Accordingly, there 
are two dimensionless parameters that describe spatial 
and temporal relationships between wave and current pa- 
rameters respectively: L = ku r is the ratio between 
the distance passed by the fluid particle during r and 
the wavelength, and T = tt k r is the ratio between the 
correlation time of currents and wave period. The char- 
acteristic function (exp [ik ■ J(t)]} generally depends on 
the details of the currents statistics but it has universal 
behavior both at t <C T and i > t where general calcu- 
lations are possible. On the plane of the dimensionless 
parameters L, T we distinguish three regions of different 
asymptotic behavior. 

Consider first the ballistic limit when the integral © 
is determined by the times t <C t when the drift ve- 
locity does not change and J(t) w u(0, x)t. Again, 
only those waves contribute that are in a Cherenkov 
resonance with the current (whose phase velocity coin- 
cides with the local projection of the current velocity): 
/i = Tr(5(Q k — fc • it)}. In this limit, the weight u is deter- 
mined by the single-time probability distribution of the 
current velocity which we denote V(u). In particular, for 
the isotropic Gaussian V(u) oc exp(— u 2 /2uq), we get 



(i(k) = (ird/2) 1/2 {ku )- 1 exp[-(Vdn k /V2ku ) 2 } , (9) 

The ballistic approximation and © hold when k 2 u^T 2 /d 
is much larger than both unity and flkT. 

The second universal limit is that of a slow cluster- 
ing which proceeds for the time exceeding the correlation 
time of currents. At t 3> r, we use the diffusion approx- 
imation, (exp [ik ■ J{t)]) = exp[— k 2 UQTt/d], in ([5]): 



d(kuQTY 



(fc Mo r) 4 + (dftfcr) 2 



(10) 



That answer and the diffusive approximation hold when 
both k 2 UgT 2 jd and VtkT are small. Formulas (|7I10|) can 
be compared with the expression for the clustering rate 
for waves with a linear damping, A ~ Jk 2 Ek^k(^k + 
r )\)~ 1 dk (2H- We see that in this limit the diffusive mo- 
tion of fluid particles due to currents is equivalent in its 
effect to a damping of waves with 7^. = k 2 u^r/d : where 
Mqt/c? is the eddy diffusivity. 

The third asymptotic regime takes place for fast- 
oscillating waves when VL^t exceeds both unity and 
k 2 UQT 2 /d. An integral of the fast oscillating expo- 
nent with a slow function, J °° cos(Qkt)f{t) dt, decays 
as Vl^ 2n ~ 2 where 2n + 1 is the lowest order of the non- 
vanishing derivative of f(t) = (exp [ik ■ J(t)]) at t = 0. 
When all odd derivatives at zero are zero, the integral 
decays exponentially. We see that the answer depends 
on the details of the statistics of currents. 




FIG. 1: The isolines of the dimensionless clustering rate 
fi(k)Q,k given by pip , here L — ku^r, T — fifcT. 



If u(t, X(t,x)) is Gaussian and isotropic with 
(u a (0,x)uP(t,X(t,x))) = (ul/d)5 a0 exp(-\t\/r) then 

u(k)=r dscos{Ts)exp[(L 2 /d)(l-s-e- s )]. (11) 
Jo 

It gives both limits (|9I10|) and 

fi(k) = (ku ) 2 /rnid, (12) 

at large flk since the lowest non- vanishing derivative is 
/'"(0). Isolines of (fTTj) are shown in Figure 1 for arbitrary 
parameters. Remind that the whole description based on 
([5]) is valid when v -C u. 

If one interpolates between the ballistic and diffusive 
regimes (i.e. between J 2 oc t 2 and J 2 oc t) with the help 
of the function a/1 + (i/r) 2 — 1, which is smooth at t = 0, 
then the weight factor can be calculated analytically: 



fj,(k)=T ds cos(Ts) exp (L 2 jd) ( a/I + s 2 - 1 



tL 2 (L 2 \ K^y/iHW+T*) 



cT exp V d 



y jL i /d 2 +T 2 



(13) 



Here K±(x) is a Bessel function of an imaginary ar- 
gument having the following asymptotics: K\(x) = 
v / n/2xexp(-x)[l + 0(1/ x)} for x > 1 and K x (x) ~ 
1/x + 0(x\n(x)) for x < 1. We see that (HSJ) repro- 
duces (0 El in the regions L 2 /d > 1, L 2 /d > T and 
L 2 /rf < 1, T <C 1 respectively. At the fast-oscillation 
limit one gets an exponentially small contribution 



'7r(fcu ) 4 - 

2d 2 n 3 k 



,-0, k T 



(14) 



That concludes the analysis of the weight fi(k) and we 
can now turn to (0) to get the clustering rate A. 

When the wave spectrum is not very wide (with the 
width comparable to k) we get 



A ~ (kvYfi{k) = e 2 nifj,(k) . 



(15) 
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Let us now find out which wavenumbers contribute ([7]) 
when the spectrum is wide. Consider an isotropic power 
spectrum oc k b ~ d between some k m j .n and k max and 
the dispersion relations Qk = Ck a [23]. Consider first 
the ballistic regime. For (a>lA6>0)V(a<lA6<0) 
the wavenumber fc* = [bu 2 /dC 2 (a — l)] 1 /( 2a_2 ) deter- 
mines A. For (b < A a > 1) V (b < -1 A a = 1) 



the surface. When compressible component of the veloc- 
ity is small, the Lyapunov exponents are due to the cur- 
rent flow, Ai ~ A2 ~ t _1 . Then, the fractal dimension of 
the density distribution can be expressed by the Kaplan- 
Yorke formula 1+Ai/|A 2 | = 2-A/|A 2 | » 2-Ar. The frac- 
tal part reaches maximum in the ballistic regime when 
Clk — kuQ, then At ~ e 2 VL k T = e 2 k£ grows with £ = uqt 



the clustering rate is determined by k m i n , while for and reaches order unity when k£ ~ e 2 . Therefore, the 



(b > A a < 1) V (b > -1 A a = 1) by fc„ 



Let 



us give physical examples using Kolmogorov spectra of 
waves. For capillary waves on a deep water, fi^ cx k 3 ^ 2 
and Ek cx fc -11 / 4 , and A is determined by k m i n i.e. by 
longest waves in the wave turbulent spectrum (assum- 
ing the ballistic approximation is valid for them). For 
gravity waves on a deep water, f2fc cx k 1 / 2 , and for both 
Kolmogorov solutions, Eu oc /c~ 20 / 6 and E\. cx k~ 7 / 2 , the 
clustering rate is determined by waves around fc* . For dif- 
fusive regime, the clustering rate is determined by k max 
if b > max[2a — 4, 0] and by k min if b < max[2et — 4, 0]. 

Estimates (|9I10I12I14|) show that \/(e 2 n k ) ~ n k fj,(k) 
is a dimensionless function which has a maximum of order 
unity either in the ballistic regime where the phase veloc- 
ity of wave is comparable to the current velocity or in the 
diffusive regime where the eddy diffusivity UqT is compa- 
rable to fifcfc" 2 (in the third asymptotic regime A/ (e 2 £lfc) 
is always small). In those cases, A/f^ ~ e 2 , i.e. the 
degree of clustering during a period is the squared wave 
nonlinearity (typically e is between 0.1 and 0.01). Such 
clustering is pretty fast (minutes for meter-sized gravity 
waves and a week for fifty-kilometer-sized inertio-gravity 
waves). Therefore, it is likely that the interplay between 
waves and currents is the source of inhomogeneities of 
floater distribution in many environmental situations. 

Clustering leads to fractal distribution of floaters over 



distribution is most fractal when waves are weak and 
short while currents are long and strong. For example, 
meter-sized gravity waves on water surface will produce 
most inhomogeneous distribution of floaters when there 
are currents with velocities in meters per second and 
scales in hundreds of meters. We see that the current- 
to-wave ratio of scales, k£, compensates for a small wave 
nonlinearity, e 2 , so that even weak waves with the help 
of solenoidal currents can produce very inhomogeneous 
fractal distribution of matter. 

As a final remark, note that apart from fluid mechan- 
ics, one can think about the evolution of a dynamical 
system as a flow in the phase space and treat density 
as a measure. Solenoidal (incompressible) flow corre- 
sponds to Hamiltonian dynamics and to a constant (equi- 
librium) measure. Compressibility corresponds to pump- 
ing and damping i.e. to non-equilibrium. Indeed, the no- 
tion of singular (fractal) measures first appeared in non- 



applied in fluid mechanics j^l^MlllT 12|. Therefore, the 



apt 

equilibrium statistical physics [1 7|. | 1S] . 19j and then was 



formulas ([6 HT4j) also describe the entropy production rate 
in dynamical systems under the action of perturbations 
periodic in space and in time. 
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